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Abstract. We study a class of mass transport models where mass is transported 
Q ■ in a preferred direction around a one-dimensional periodic lattice and is globally 

conserved. The model encompasses both discrete and continuous masses and parallel 
and random sequential dynamics and includes models such as the Zero-range process 
and Asymmetric random average process as special cases. We derive a necessary and 
sufficient condition for the steady state to factorise, which takes a rather simple form. 
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Mass transport models form a general class of lattice models defined by dynamics 
in which mass is transferred (without loss) stochastically from site to site. They 
have attracted much recent attention, especially in connection with "condensation 
transitions" [U 121 El IU El H3 • Examples include the Zero- Range Process (ZRP) pQ and 
Asymmetric Random Average Process (ARAP) IB] , which have been used to model 
such diverse situations as traffic flow, clustering of buses [3j, phase separation dynamics 
and force propagation through granular media [10]. In general, it is difficult to 
determine the steady state distribution of such models. Thus, it is remarkable that, not 
only are the steady states of many models found, they often share a very convenient 
property, namely, a factorised steady state (also referred to as a product measure). Of 
course, such a property greatly facilitates the analysis of interesting behaviour, e.g., 
condensation. 

In this letter we determine the requirement for a factorised steady state in a 
very broad class of mass transport models. The form of this necessary and sufficient 
condition, stated in (|TH|) . turns out to be appealingly simple. Encompassing both 
random sequential and parallel dynamics, this class includes both the ZRP and ARAP. 
We discuss the salient features of the approach leading to (fTKJ) and recover some 
previously known cases. 

We consider a one-dimensional lattice of L sites with periodic boundary conditions 

(site L + 1— site 1): associated with each site is a mass rrii, i = 1 . . . L. The total mass 

is given by M = Yli=i m i- We shall most generally consider rrii as continuous variables. 

The dynamics is as follows: from time t to t + 1, at each site %, mass //j drawn from a 

distribution </?(//j|mj) 'chips off' the mass rrii, and moves to site i + 1. Thus the master 

equation for the weights (unnormalised probabilities) F t (rn) is 

£i /"oo r m 'i 

F t+1 (m) = T [ / dml / d/^ <f{jii\m'^ TT 8{m~m' j J rHj-Hj- 1 ) F t (mf_) , (1) 
i=1 Jo Jo j=1 

where m = {mi, m 2 , ■ ■ ■ , m L }. Note that this dynamics conserves the total mass, M, so 
that F t (m) may be considered as a function of only L — 1 variables. The integral of the 
weights, subject to the constraint of globally conserved mass, 

Z{M,L) = ]JI drriiSlM-J^mjFtim), (2) 

should be finite and serves as a "partition function," so that F/Z is a probability density 
(or "canonical distribution"). 

In the t — > oo limit, F t (rn) approaches a ^-independent function, which we denote 
simply by F(m) and refer to as the steady state. Defining the Laplace transform 
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F{m) , (3) 



and transforming ([!]). we find 
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JT / drn^ / d^ p(^|mj)e-* i ( m *- w+w - 1 ) 
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F(m') . (4) 
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We now assume that the steady state weight factorises 
F(m) = Hf(m i ) 

which implies 

/•oo 

G(s) = Y\g{si) where g(s) = / dmf{m)i 
i J o 

Then (@J becomes 

n^)=n 
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Changing variables to o = m — /i (the mass remaining after the move), we write 

f(m)<p(ji\m) = V(ji,<r). (8) 

Note that no assumption on the form of f(m) or yp(/j,\m) is implied at this point. With 
this notation (J7J) becomes 

r /*00 /*00 

i i L-/0 Jo 

A necessary and sufficient condition for the solution of (|5j). is 

/*oo /»oo 

/ dm da i V{fx i ,a i )e- s ^ ai +^^£{s i )k{s i+1 ) 
Jo Jo 

where the two functions, k and £, must satisfy 

k(s)e(s)=g(s). 

That ()10j) is necessary and sufficient may be seen by taking the logarithm of Q then 
taking derivatives with respect to s, and Si+%. 

Condition (fTT|) implies via the convolution theorem that 



(9) 

(10) 
(11) 



/ (m) = [v * w] (m) = / d/j, v (fi) w (m — //) 



where 



fcOrt 



cfyie SM w(yu) 



£(s) = / dae- sa w(a). 
Jo 



Then, equations (|TU|) and (|T3|) imply 
P (/i, a) = v (/i) u> (a) . 

Finally we obtain from (jHJ) and |T2|) 

v{ij)w{m-n) 
cp(jj.\m) = 



(12) 
(13) 
(14) 
(15) 



[t> * tu] (m) 

Let us emphasize that the condition for a factorised stationary distribution for the 
whole lattice precisely reduces to the condition that ip(n\m) has the form ()15|). Thus, 
equation (JT5|l is the central result of this paper: for chipping rules of the form ([T5]). 
one has a factorised steady state © with weights given by (|T2|l. Let us comment on 
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several important points. Equations (|11H15|) allow us to define "equivalence classes" 
of chipping distributions — those leading to the same stationary state — by dividing k (s) 
and multiplying I (s) by any (well behaved) function of s. In particular, their roles can be 
"reversed" to form a "dual" ip, i.e., w (fi) v (m — //) / [w * v] (m). Now, we can obviously 
interpret the factors in (fTlJ) as a function for //, the mass which moves, and a function of 
a, the mass which stays. In this sense, "duality" reverses these two portions of the mass, 
without changing F(m). If we further perform a Galilean transformation (shifting the 
entire lattice by one site in a time step) and a parity transformation (i -£=> L + l — i), then 
we recover the original system. Finally, note that both ip and the steady state {F/Z) 
are invariant under shifts of lnt> and In w by a linear function (i.e., there are arbitrary 
amplitudes or exponential factors in v and w: a M , a a ). 

In addition to treating models with parallel dynamics, manifest in (JTJ), we can 
extend the approach outlined above to models with random sequential dynamics. Let 
the probability of a chipping event in a time step ex dt so that, to leading order in dt, 
at most one chipping event over the whole lattice occurs per update. Furthermore, we 
can let the duration of a time step be dt and take dt —>■ to obtain a continuous time 
limit where chipping events occur with rates per unit time. Thus, we write 

v (/x) = 5(fi) + x(fi)dt, (16) 

where <5(/i) is the Dirac delta function. Then (|T2j) and ([To]) yield 

f(m) = w(m) + dt[x * w](m) and 

ip(u\m) = — ; — r — — — - \5(u)w(m) + dtx(u)w(m — u)} 

rv ^' ' w(m) + dt[x*w]{m) 
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Taking dt — > we obtain the continuous time limit where mass jjl moves from a site 
with mass m with rate x(fj,)w(m — n)/w(m) and f(m) = w(m). 

Let us illustrate how this approach unifies two seemingly unrelated models - ARAP 
and ZRP. First we consider the ARAP [3IBIEJEJ, a model in which each site contains 
a continuous amount of mass and at each time step a random fraction of the mass 
moves to the next site to the right. Its precise definition lies in i[)(r\m) = (p(fi\m)m, the 
distribution for r, the fraction of mass that moves to the neighbouring site. A known 
family of distributions where one has a factorised steady state is ip(r\m) = (n — l)r n ~ 2 
[TIH ITT] which becomes 

^H = (n-l)4rr- (17) 

m n 1 

In our approach, the results are particularly simple: 

v(ji) =/i"- 2 , w(a) = l, (18) 

f( m ) = m"" 1 / (n - 1) . (19) 

Note that, to relate this / (m) to relevant quantities in the literature (e.g., [UJ), the 
single site mass distribution, defined as the full distribution integrated over the rest 
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of the mass variables, is p(m) = f (m) Z (M — m,L — 1) jZ (M,L). In this case, 
Z (M, L) = M nL ~ 1 [r(n — 1)] /r(nL), so that our p (to) reduces, e.g., to equation (37) 
of ^I] in the thermodynamic limit. 

Another well-known case is the Zero- Range Process, reviewed in pQ. A focus of 
major interest (for recent developments see for example fHJ O E3 E3), it is a mass 
transport model where m^ takes integervalues and a unit mass moves from site i to site 
i + 1 with probability u(mi). Within our approach, this model appears as a very special 
case, with Dirac delta distributions for both v and w. Since the moved mass can take 
only two values while the one remaining can be of any integer, the most general forms 
are 

oo 

v(ii) = <5(/i) + a<5(/i — 1), w(o~) = y^Wkd~(a — k) , (20) 

A:=0 

where a and w^ are arbitrary weights. As overall amplitudes are irrelevant, we have 
chosen the coefficient of 5 (fi) to be unity and will set Wq = 1. From / = v * w, we see 
that 

f(m) = w(m) + aw(m — 1) (21) 

oo 

= 5{m) + y^ [aw m _i + w m ] 5(m — k) . (22) 

fc=i 
With a little care, we obtain 

w m 5(fi) + aw m -. 1 6(jj, - 1) 

cp(jJL\m) = . (23) 

w m + aty m _i 

The coefficient of S(fi — 1) is precisely the chipping probability, denoted by u (to) above. 
From here, we easily find the w m in terms of the u: 



. 1 - u(n) 

w m = a 

■ ■ u{n) 

n=l 



nL — mill , 

, • 24 

u(n) 

ra=l v ' 

Substituting this expression into (|2*2*jl yields for the weights, 

oo 

/H = J]/ fc (J(m-fe) (25) 



fc=0 



where 



/^r-^n 1 ^ &» ^i, (26) 

1 — m (/c) - LJ - u(n 

v ' n=l v ' 

= 1 for k = . (27) 

This result was previously obtained by a more complicated approach ^7j. Note that 
the factors a k will drop out when we consider the probability density itself: F(m)/Z. 
We close this paragraph by noting the case with random sequential dynamics, which is 
obtained by letting a = adt and u(m) = x(m)dt where dt — > yielding 

k 
n=l v ' 
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Finally, the results presented here may be generalised to the case of heterogeneous 

mass transfer where (pi(fj,\m) depends on the site i. A necessary and sufficient condition 

for a factorised steady state is that 

, , . v (u)wi(m — a) ,_. 

^(fxm) = 7 l \ fl 29 

[v * Wi\ (m) 

where v and Wi are arbitrary functions but v must be the same for each site. The weight 

functions are given by 

fi(m) = [v*w l ] (m) . (30) 

To conclude, we have determined the condition for the steady state in a general 
class of mass transport models to factorise. This class encompasses both continuous 
and discrete mass, as well as parallel and random sequential dynamics. Not only does 
this approach provide a unified perspective of all previously known models, it opens 
avenues to construct new models with this property (e.g., binomial chipping process 
and generalized Zero- Range Processes). In addition, we believe this approach would 
facilitate a deeper understanding of the existence and nature of condensates and possibly 
reveal novel forms of phase transitions. Implications of the gauge-like transformations 
should also be explored. Further work is in progress and will be published elsewhere. 

Acknowledgments 

We thank Max Planck Institute, Dresden where this work was initiated for hospitality. 
One of us (RKPZ) thanks H.W. Diehl for his hospitality at the University of Duisburg- 
Essen, where some of this research was carried out, and acknowledges the support 
from the Alexander von Humboldt Foundation and the US National Science Foundation 
through DMR-0088451. 

[1] M.R. Evans, Braz. J. Phys. 30, 42 (2000) 

[2] M. R. Evans, Europhys. Lett. 36, 13 (1996) 

[3] P. Bialas, Z. Burda, and D. Johnston, Nucl. Phys. B 493, 505 (1997). 

[4] O.J. O'Loan, M.R. Evans and M.E. Cates, Phys. Rev. E 58, 1404 (1998) 

[5] S.N. Majumdar, S. Krishnamurthy and M. Barma, Phys. Rev. Lett. 81, 3691 (1998); J. Stat. Phys. 
99 (2000) 1. 

[6] E. Levine, D. Mukamel, G. Ziv (2004) Condensation transition in zero-range processes with 
diffusion cond-mat/0402520 

[7] J. Krug and J. Garcia, J. Stat. Phys., 99 31 (2000) 

[8] R. Rajesh and S. N. Majumdar, J. Stat. Phys., 99 (2000) 943; Phys. Rev. E 64 (2001) 036103. 

[9] Y. Kafri, E. Levine, D. Mukamel, G.M. Schutz and J. Torok, Phys. Rev. Lett. 89, 035702 (2002) 

[10] S.N. Coppersmith, C.-h. Liu, S. Majumdar, O. Narayan, T.A. Witten Phys. Rev. E., 53, 4673 

(1996) 

[11] F. Zielen and A. Schadschneider, J. Stat. Phys 106, 173 (2002) 

[12] F. Zielen and A. Schadschneider, J. Phys. A 36, 3709 (2003) 

[13] M. R. Evans and T. Hanney, J. Phys. A 36, L441 (2003) 

[14] S. GroBkinsky, G.M Schutz and H. Spohn, J. Stat. Phys 113, 389 (2003) 

[15] C. Godreche, J. Phys. A: Math. Gen. 36, 6313 (2003) 

[16] A. G. Angel, M. R. Evans and D. Mukamel, J. Stat. Mech.:Theor. Exp. P04001, (2004) 

[17] M. R. Evans, J. Phys. A: Math. Gen. 30 5669 (1997) 



